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Experimental studies on ion-water clusters have provided insights into the microscopic aspects 
of hydration phenomena. One common view is that extending those experimental studies to larger 
cluster sizes would give the single ion absolute hydration free energies not obtainable by classical 
thermodynamic methods. This issue is reanalyzed in the context of recent computations and molec- 
ular theories on ion hydration, particularly considering the hydration of H + , Li + , Na + , and HO~ 
ions and thence the hydration of neutral ion pairs. The hydration free energies of neutral pairs 
computed here are in good agreement with experimental results, whereas the calculated absolute 
hydration free energies, and the excess chemical potentials, deviate consistently from some recently 
tabulated hydration free energies based on ion-water cluster data. We show how the single ion 
absolute hydration free energies are not separated from the potential of the phase in recent analyses 
of ion- water cluster data, even in the limit of large cluster sizes. We conclude that naive calculations 
on ion-water clusters ought to agree with results obtained from experimental studies of ion-water 
clusters because both values include the contribution, somewhat extraneous to the local environment 
of the ion, from the potential of the phase. 
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I. INTRODUCTION 

The interactions between the ions and inner shell wa- 
ter molecules are much stronger than thermal energies 
and typically display chemical complexities even in the 
simplest cases 0, Q . Computational chemistry and the 
theories of solutions are approaching the stage that calcu- 
lations on the free energies of ions in water must be taken 
seriously. Despite this progress, there is still a non-trivial 
level of confusion about what is to be learned from exper- 
imental comparisons, particularly with comparisons de- 
veloped from current experimental studies of ion-water 
clusters. This paper attempts to provide guidance on 
these issues by presenting a thermodynamic reconsidera- 
tion of conclusions drawn from ion-water cluster experi- 
ments, and results of recent molecular calculations. 

Computational results are presented for aqueous ions 
H + , Li + , Na + , and HO~ which participate in a myriad 
of biological processes Q. These are common ions in 
aqueous phase chemistry, and this is especially true for 
H + and HO~. These species also play a significant role 
in current problems such as the speciation of beryllium 
metal in the environment or in the lungs. Under- 
standing beryllium speciation, and the development of 
beryllium toxicity, is of substantial technological inter- 
est, and would benefit from a molecular understanding 
of ion hydration phenomena. 

Absolute hydration free energies, or equivalently sin- 
gle ion activities, that we wish to clarify are not mea- 
sureable by purely thermodynamic means 0, |(|; cor- 
rect single ion activities will contain fundamental extra- 
thermodynamical information and hence their useful- 
ness in thermodynamical analysis might be questioned. 
Extra-thermodynamical, single ion activities, however, 
are not unknowable (see |7| and Furthermore, that 
knowledge would clarify the molecular understanding of 
ion hydration, and permit a compact, tabulated consen- 



sus of molecularly valid, experimental thermodynamic in- 
formation. 

The existing tables, for example [HEJEIi should nat- 
urally be consistent in describing hydration free energies 
for neutral combinations of ions, since those combina- 
tions are thermodynamically measurable. But the align- 
ment of those tables depends on extra-thermodynamical 
information. An accurate determination of the properly 
defined single ion activity for any ion represented in those 
tables would be sufficient to align those tables. It has 
been pointed out recently 0] that in favorable cases the 
inaccuracies of computed single ion hydration free ener- 
gies could probably be made less than the misalignments 
in current tables of single ion hydration free energies. 
The present paper pursues this possibility further. 

Because simulated extra-thermodynamical informa- 
tion is available from molecular simulations, it is possible 
to compute solely the work of coupling the ion with the 
solvent in its locality, i.e., the absolute hydration free 
energies of the ion separated from the contribution from 
the potential of the phase (section II). Simulation results, 
however, critically depend on parameterized molecular 
models and on the treatment of electrostatic interactions. 
Moreover, for many ions, particularly H + and HO~, the 
interactions with water are not easily described by con- 
venient force-fields. This problem is circumvented by the 
quasi-chemical theory of solutions, within which the ex- 
cess chemical potential of the solute is partitioned into an 
inner-sphere contribution, accounting for the chemistry 
in the inner shell, and an outer-sphere contribution. This 
latter contribution describes the interaction of the inner- 
sphere quasi-component with the rest of the fluid and 
can usually be well described by simple force-field mod- 
els or dielectric continuum models. This quasi-chemical 
approach has been used before to treat ions, such as, H + 
PJ, Li+Q, Na+0, and HO" [lj|. Chemical reactions 
in solution, such as speciation of Fe 3+ and Be 2+ ^(| , 
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are also well described. 

In this article we bring together our earlier observa- 
tions on the hydration of monovalent ions in conjunction 
with classical force-field based estimates of outer-sphere 
contributions. The results so obtained for neutral ion 
pairs are in good accord with experiments, but our sin- 
gle ion values differ consistently from current tabulations 
based on cluster-ion hydration experiments. 

Our plan for this report is as follows: First, we define 
our problem and establish the notation we will use. We 
then discuss the experimental studies leading to current 
estimates of absolute hydration free energies of the ions. 
After that, we present the quasi-chemical theory and our 
computational results. Finally, we return to analyze the 
experimental cluster-ion hydration results in this context. 
In the last section we identify conclusions of this work. 



II. ABSOLUTE HYDRATION FREE ENERGIES 
OF IONS IN WATER 

We will consider the chemical potentials of ionic solutes 
in liquid water and we cast these chemical potentials in 
the form 

Mm+M = e<i>{w) 

+ RT In [p M+ HA M+ 3 ] 

+ Mm+H ■ (!) 

This treats an ionic solute M + suggested to be a metal ion 
of charge e, but the notation here will be extended nat- 
urally to anions X - . The qualifier [w] in Eq.^indicates 
the macroscopic phase to which this chemical potential, 
or another quantity there, applies. Thus, Pm+( w ) is the 
number density of M + in a liquid water phase. (Though 
we will be particularly interested in the case of liquid 
water as a solvent, other cases must be permitted also 
in the typical thermodynamic analysis.) The quantity 
4>(w) is the electrostatic potential of the w phase. This 
may be introduced on the basis of the principle that a 
mean electric field must be zero through interior regions 
of a macroscopic conductor at equilibrium, and we will 
use the language that 4>(w) is the potential of the phase. 
Only differences in these potentials between phases will 
be involved, however, in our analysis. The quantity A M + 
is the thermal deBroglie wavelength for the species M + , 
a known function of the internal characteristics of an M + 
ion, and of the temperature T, but with no dependences 
of further thermodynamic relevence. (The standard state 
adopted by Eq. ^ is sometimes referred to as the Ben- 
Nairn standard state.) 

Finally, ~p M + (w) is the object of the present study. It 
depends on temperature, pressure, and composition of 
the system. In the limit of vanishing solute concentra- 
tion we will call it the absolute hydration free energy of 
M + . We are principally interested in standard temper- 
ature and pressure conditions. As Eq. ^ is formulated, 
Mm+ ( w ) depends on the energetic interactions of solution 



constituents, i.e., the solvent and dissolved species, with 
M + ions, and would vanish if those interactions were to 
vanish. That is also the case for the contribution e<f>{w). 

The separation imposed by Eq.^Jis clear on the basis of 
our mechanical understanding of electrostatics. But from 
a thcrmodynamical view it is a nonoperational separation 

m 

As a gauge of molecular structure and energetics of the 
solution in the neighborhood of an M + ion, ~p M + (w) is 
worthy of interest. In contrast, <f>(w) reflects distribution 
of charge on the surface of the system, and is of less 
intrinsic interest because it is not sensitive exclusively 
to the local condition of an ion of interest. The goal is 
to separate /Im+ ( w ) horn thermodynamic combinations 
such as 7Im+( w ) + ~Px- ( w )> to which the potentials of the 
phases don't contribute, and to separate ~p M +(w) from 
e<p(w). 

Our notation here is deficient in one respect. For equi- 
librium of two conducting phases containing ionic solutes 
at low concentration, the potential difference A<p = ^(7) 
- 4>(rf) across the phase boundary is not simply a prop- 
ertyof the pair of solvents in contact with one another 
[a (Note that many materials can be considered 

conductors for sufficiently long times.) From a formal 
point of view of macroscopic thermodynamics, this fol- 
lows from the requirement that the bulk compositions of 
a macroscopic conductor be electrically neutral. More 
physically, the ions make a Donnan-like contribution to 
the difference between the electrostatic potentials of the 
phases to achieve neutral bulk compositions. This poten- 
tial difference does not progressively vanish as the con- 
centration of the ionic solutes decreases, so it is not op- 
erationally avoidable by confining attention to low elec- 
trolyte concentrations. The difference A<f>, whatever its 
source, does not contribute to the thermodynamic combi- 
nations /I M + (w) + 7*x- ( w )- Therefore an accurate evalu- 
ation of 7*m+ ( w ) on an y basis would avoid consideration 
of A<fi altogether. 

Thermodynamic considerations will involve differences 
in these free energies and differences A(f> in the electro- 
static potentials. If a dilute gas is considered to be the 
coexisting conducting phase, it is instinctive to assign the 
potential of that phase to be zero. No other considera- 
tions need be affected by this choice. It is then convenient 
to take <j>(w), for example, to be the potential relative to 
a dilute gas phase. 

III. ANALYSES BASED UPON ION- WATER 
CLUSTER DATA 

Discussions of the difficulty of obtaining the absolute 
hydration free energy ( w ) often begin with the ob- 
servation that thermodynamic processes for bulk phases 
involve manipulation of neutral combinations of material. 
It is natural, therefore, to consider adding a single ion to 
a sub-macroscopic amount of solvent, for example 

M+ + (H 2 0)„^M(H 2 0)„+ (2) 
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for n not too large. This leads to the consideration of 
ion- water cluster experiments as a potential source of in- 
formation on ~p M +(w). Coe 0] helpfully reviews that 
work. These experimental methods permit determina- 
tion of a standard contribution to the chemical potential 
of the cluster on the right side of Eq.[21in a dilute gas, for 
cluster sizes typically including n=4, 5, 6. It is helpful 
to keep in mind that these quantities can be targets of 
molecular simulation j^jj, EI as well. 

We collect several relevant points about these quan- 
tities. A first point is that the standard chemical po- 
tential of such a cluster includes contributions from the 
kinetic motion of the ion over the interior of the clus- 
ter. For mesoscopic clusters, this contribution may be 
awkward because the volume of the cluster accessible 
to the ion can be ambiguous. An extreme example of 
this ambiguity is the p henomenon , no w well appreciated 

m m m eeihiiiiiiiiiim mm m, that some 

ions might be localized on the surfaces of small water 
clusters. Fortunately, the cluster sizes studied in ion- 
water cluster experiments considered here are so small 
that even the surface vs. interior question doesn't ap- 
pear to be a significant problem. 

A second basic point is that neither the experiment 
nor the typical simulation calculations further separate 
A/it M+ (u>) = ecf) + /I M + (w) into an absolute hydration 
free energy and a contribution from the potential of the 
phase. This is natural because both of these contribu- 
tions derive from intermolecular interactions, and in the 
case of simulations it would take a special separate calcu- 
lation to evaluate the potential of the phase separately. 

Our final basic point about the separation Eq. ^ is 
that the difficulties in separating the various contribu- 
tions to the interaction contribution to the chemical po- 
tential Eq. ^ from physical data are not trivially simpler 
for a particular n. Though the contribution e(j) is ex- 
pected to depend on n, it will typically vanish neither for 
small nor for large values of n. 

Returning to the experimental analysis, we note that 
the conventional hydration free energies of ions are ref- 
erenced to the case for H + (aq). So the conventional hy- 
dration free energy for the M + ion is ~p M + (w) — ~p H+ (w) 
and similarly that for X - is /I x _ (w) +7%+ i w )- Following 
Klots [34j . we then form the following combination and 
obtain: 

A M M + ( w ) - A ^x+ ( w ) = 2(e0(io) + 7Z H+ (id)) 
+ (7%+M -Mh+H) 
- (Px-N +Mh+H) • ( 3 ) 

The last two contributions in parenthesis on the right 
are conventional single ion hydration free energies and 
can be obtained from tables. Klots gives a plausible 
procedure for obtaining the quantities on the left from 
the ion-cluster experimental results. If the quantities 
on the left of Eq. [3| were obtainable to sufficient accu- 
racy from cluster studies (Eq. [2) , then Eq. could be 
solved for e<p(w) +Jl- [i+ (w) = A^ h+ (uj), corresponding 



to the same standard state as for H + in compiling the 
conventional ion hydration free energies. The evidence 
presented 0, 0] suggests that this procedure is suc- 
cessful to an interesting degree, and remarkably so for 
small values of n. In view of those results and the analy- 
sis from Eq. [21 we conclude that the values otained have 
not separated the absolute hydration free energy ~p M+ (w) 
from the potential of the phase in A/i H+ (w). 

In concluding this section, we return to consider the 
Donnan-like contribution to A0 that was noted above 
[H 0| . This electrostatic contribution has its source in 
sorting of ions of different charges within the ion corre- 
lation length of an interface. Though the source of this 
contribution is the ions distributed in interfacial regions, 
this contribution typically achieves a non-zero value in 
the limit of vanishing ionic contributions. This contri- 
bution is present in the macroscopic description Eq. ^ 
But because the ion-water cluster experiments consid- 
ered here treat clusters with only one ion that Donnan- 
like ionic effect isn't present in those experimental data. 

IV. QUASI-CHEMICAL THEORY 

In the quasi-chemical theory |35| . the region around 
the solute of interest is partitioned into inner and outer 
shell domains. A variational check of this partition ex- 
ists, if desired U|. Here we treat the inner shell, where 
chemical effects are important, quantum mechanically. 
The outer shell contributions have been assessed using 
both a dielectric continuum model and classical molecu- 
lar dynamics simulations. 

The inner shell reactions pertinent to the hydration of 
X ± ( = H+, Li+, Na+, HCT) are 

X ± +nH 2 O^X - [H 2 0]„ ± 

Based on earlier work, the inner shell hydration num- 
bers are n = (2,4,4,3) for H+ Q2, Li+ 0, Na+0, and 
HO ~ 0| j respectively. The free energy change for these 
reactions were calculated using the Gaussian suite of pro- 
grams 36]. The X ■ [H 2 0] n ± clusters (and H 2 0) were 
geometry optimized in t he g as phase using the B3LYP 
hybrid density functional |37| and the 6-31+G(d,p) basis 
set. Frequency calculations confirmed a true minimum, 
and the zero point energies were computed at the same 
level of theory. Single point energies were calculated us- 
ing the 6-311+G(2d,p) basis set. 

Dielectric continuum model: For estimating the 
outer shell electrostatic contribution, the ChelpG method 
[38| was used to obtain partial atomic charges. The non- 
electrostatic contributions are expected to make negligi- 
ble contributions to the hydration free energy and are not 
considered further. Then with the radii set developed by 
Stefanovich et al.^, surface tessera were generated [iof . 
and the hydration free energies of the clusters were calcu- 
lated using a dielectric continuum model . With this 
information and the binding free energies for the chem- 
ical reactions, a primitive quasi-chemical approximation 
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to the excess chemical potential of X 1 * 1 (aq) in water is: 



PP. 



X±(iu) 



In (K n pa a o n 



(4) 



where K n = K n ^ exp 



-P (mx(h 2 o)„±H - ™/% 2 oH 
K^ ' is the equilibrium constant for the reaction in an 
ideal gas state, with n of Eq. 01 the hydration number of 
the most stable inner shell cluster, and (3 = 1/keT. The 
density factor ph 2 o appearing in eq.Q] reflects the actual 
density of liquid water and its effect is accounted for by 
including a replacement contribution of — nkeT ln(1354). 
A detailed statement on standard states and this replace- 
ment contribution can be found in Grabowski et al. |l2j| . 
Relevant energies are collected in TABLE |T| For the 



TABLE I: Electronic energy (a.u.), thermal corrections (a.u.) 
to the free energy, and electrostatic contributions to the 
excess chemical potential (kcal/mole) using dielectric con- 
tinuum (DC) approximation and molecular dynamics (with 
TIP3P and SPC/E water models). The statistical uncer- 
tainties in the molecular dynamics values are of the order 
of 1.5 kcal/mole (Fig.0. Lhe partial charges are obtained at 
B3LYP/6-311+G(2d,p). 





E 


G corr 


DC TIP3P SPC/E 


HaO 


-76.45951 


0.00298 


-7.7 


-6.6 


-6.5 


H+ 





-0.01 








H[H 2 0] 2 + 


-153.24860 


0.03253 


-77.5 


-66.0 


-64.5 


HO" 


-75.82779 


-0.00771 








HO[H 2 0] 3 - 


-305.32036 


0.04705 


-72.7 


-85.2 


-87.6 


Li+ 


-7.28492 


-0.012748 








Li[H 2 0] 4 + 


-313.29314 


0.05758 


-64.0 


-54.2 


-51.4 


Na+ 


-162.08757 


-0.014429 








Na[H 2 0] 4 + 


-468.05512 


0.05095 


-62.0 


-51.5 


-49.8 



dielectric continuum calculations, an observation volume 
of radius 1.9 A and 2.3 A was centered on the Li + and 
Na + ions, respectively. For H + and HO + these radii 
were 1.172 A and 1.56 A respectively. The calculated 
values change only modestly to increases in these values. 

Classical molecular dynamics: The electrostatic 
contribution to the hydration (excess) free energy is given 
by 



Jo 



(5) 



where A is the coupling parameter which switches the so- 
lute charge from to q. (ip)\ is the ensemble-averaged 
potential on the ion at a particular charge state, A. 
Gauss-Legendre quadratures provides a facile route 
to estimate this integral. In particular, treating (ip)\ as a 
third degree polynomial, the following two-point formula 



Wa_) 



(6) 



where A± = q(l/2±l/\/l2), is exact to 4 th order in per- 
turbation theory |42j| . Trial calculations using a purely 
linear response approximation (2 nd order in perturba- 
tion theory) deviate modestly from the above approxi- 
mation, and hence for greater realism we employed the 
above form. For clarity the above expressions pertain to 
an atomic ion, but are easily generalized to complex ions 
with distributed charges. 

Since the solute is largely buried within the hrst shell 
water molecules, there is no need for developing accurate 
parameters for those. This is true for the cations. For 
HO~, the hydroxyl hydrogen is largely free ^4|. But the 
charge on this hydrogen is small, being about half the 
charge of the hydrogen atoms in the classical water mod- 
els [lj|- Based on these observations, the van der Waals 
parameters for Li + and Na + were obtained from |43| and 
[44j , respectively. The H and O atoms in the cluster were 
assigned the same van der Waals parameters as those in 
classical water model, of which both TIP3P[45j |4(j and 
SPC/E 47] were used. The charge distribution on these 
clusters were the ChelpG charges. The clusters were sim- 
ulated in simulation cells of different sizes to estimate ion 
finite-size effects, and included 32, 64, 216 and 306 water 
molecules in addition to the cluster of interest. 

The simulations were all performed with the NAMD 
|48j code which uses the particle mesh Ewald method to 
treat long range electrostatic interactions. For each A, 
the simulation consisted of 200 ps of equilibration with 
velocity scaling every 250 fs, followed by another 200 ps 
of equilibration with velocity scaling every 1 fs. The wa- 
ter geometry was constrained using SHAKE. The clus- 
ter was held fixed. For the 32 and 64 water molecule 
boxes, a production run of 375 ps was conducted. For 
the 216 and 306 water molecule boxes, the production 
run lasted 150 ps. Frames were stored every 5 fs. Then 
using programs developed in-house, the Ewald potential 
at the solute charge sites were calculated, and statistical 
uncertainties were estimated by block-averaging in block 
sizes of 7.5 ps, corresponding roughly to the dielectric 
relaxation of pure water |49j. Then the free energy was 
assembled. 

Electrostatic finite size corrections were applied as de- 
scribed in Hummer et al. |2l|. As in Hummer et al. [2lj| . 
the real-space screening factor was rj = 5.6/ L, where L 
is the box length. A cutoff of k 2 < 38(2tt/L) 2 was ap- 
plied in Fourier space, resulting in 2 x 510 k vectors. A 
cutoff of L/2 was applied for the real-space electrostatic 
interactions. 

To check the accuracy of our code, we simulated some 
earlier results for ion hydration, including the charg- 
ing free energy of imidazolium 12 lj . Hummer et al. 
[2l| studied imidazolium charging using Monte Carlo 
techniques in boxes ranging from 16 water molecules 
to 512 water molecules. In the 16, 64, and 512 water 
molecule boxes, they find that the average interaction en- 
ergy at full charge is — 119.7 ± 0.24 kcal/mole, —113.3 ± 
0.36 kcal/mole, and —110.3 ±0.4 kcal/mole, respectively. 
We find -118.9 ± 1.6 kcal/mole, -113.6 ± 1.3 kcal/mole, 
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and —111.3 ± 2.2 kcal/mole, respectively. Given the dif- 
ferent methodologies this agreement is good. 

Figure ^ shows that for the clusters under considera- 
tion, ion finite size effects are quite modest. The values in 
table [I] using the TIP3P water model are obtained from 
fig. n] The least square line fits are weighted by the sta- 
tistical error. Since the ion finite size effect is small, the 
TIP3P and SPC/E values (TABLED were all obtained 
with the box of 306 water molecules. 
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FIG. 1: Electrostatic contribution to the excess chemical po- 
tential of the clusters using the TIP3P water model. L is the 
box size and from left to right, the boxes contain 306, 216, 
64, and 32 water molecules, respectively. The box volume 
is adjusted to give a density of 1 gm/cc including the water 
molecules part of the cluster. The ions are given a partial 
specific volume of zero (0), but this choice is more for con- 
venience and a more appropriate negative value would not 
change the results substantially. 



V. RESULTS 

The hydration free energies obtained with the three dif- 
ferent calcualtions of the outer-sphere contributions are 
given in TABLE[H] The agreement between our absolute 
hydration free energies and Coe and coworkers' 0, 
values is poor for the molecular dynamics results, and 
fairs only somewhat better for the dielectric continuum 
results. The computed hydration free energies for neutral 
ion pairs, however, are in markedly improved agreement 
with the experimentally suggested values table ITTT1 par- 
ticularly for the molecular dynamics vales which are in 
error by about 2-3 kcal/mole compared to the dielcetric 
continuum predictions which err by 8-11 kcal/mole. 

pK w of water: A further test of the pair hydration 
free energy for H + and HO~ is to compute the pK w , a 
quantity of immense interest in solution chemistry, espe- 
cially protein biochemistry. 

Earlier, Haymet and coworkers [5ll, Guissani and 
coworkers [5^] and Tawa and Pratt l53l considered the 
ionization of water. Tawa and Pratt 1531 went further, 



TABLE II: Excess chemical potentials (kcal/mole) of the so- 
lutes obtained with outer-sphere contributions based on the 
dielectric continuum (DC) and classical molecular dynamics 
approaches. The experimental values suggested by Coe et 
al. |5(} pertaining to the transfer of the solute from 1 atm 
(ideal gas) to 1 M (ideally diluted solute) is included after 
adjustment for 1 M (ideal gas) to 1 M (ideally diluted solute) 
transfer. 



DC TIP3P SPC/E Expt 



H+ 


-254.6 


-245.3 


-244.0 


-265.9 


Li+ 


-120.5 


-115.1 


-112.7 


-128.4 


Na+ 


-96.1 


-90.0 


-88.7 


-103.2 


HO" 


-105.3 


-121.1 


-123.8 


-104.9 



TABLE III: Solvation free energy of neutral ion pairs 
(kcal/mole). The solutes are transferred from 1 M (ideal gas) 
to 1 M (ideally diluted solute). Experimental numbers are 
from table [n] above, t For water, based on the experimen- 
tal gas phase free energy of dissociation (384.1 kcal/mole), 
the known pK w of water (15.7) and the known hydration 
free energy of water (—6.3 kcal/mole), we calculate a value 
of —366.6 kcal/mole for the HOH pair. Klots |34l quotes a 
value of —368 kcal/mole for this quantity. 



DC TIP3P SPC/E Expt 



HOH -359.9 
LiOH -225.8 
NaOH -201.4 



-366.4 
-236.2 
-211.1 



-367.8 -370.7| 
-236.5 -233.3 
-212.5 -208.1 



describing this ionization at different thermodynamic 
states. These are important works, but are not ab initio 
descriptions of dissociation. 

The reaction under consideration is the ionization of 
water (in water solvent): 



H 2 



IT 



The equilibrium constant for the reaction can be written 
as 



v qH+qHO- / / \ 
K w = WPo) 

qHOH 



(7) 



Here p =l M is the reference concentration for H + and 
HO~. /?=55 M is the reference concentration of water, 
qi = qi lg ((exp _/3AU )) is the partition function. qj lg is 
the ideal gas partition function, and ((exp~^ AU )) is the 
Widom factor. AU is the interaction energy of the solute 
with the solvent jH^]- — RTln((exp~ /3AU )) is precisely 
the excess chemical potential obtained above using the 
quasi-chemical approach. 

For the ionization reaction in gas phase, the free 
energy change is computed (using TABLE P) to 
be 383.4 kcal/mole, which compares favor ably with 
384.1 kcal/mole estimated experimentally (5^,|56j]. The 
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calculated pK w 's are 21.1, 15.6 and 14.5, respectively, 
for the outer-sphere contributions using dielectric contin- 
uum, TIP3P water, and SPC/E water. The agreement 
between the molecular dynamics results and the experi- 
mental value of 15.7 [53 is excellent. 



VI. DISCUSSION 

From a molecular standpoint, the quantity of most in- 
terest is /I x , the absolute hydration free energy. This 
absolute quantity is understood to be relative to the av- 
erage potential of the phase, or equivalently to the value 
of that potential of the phase being taken as zero. 

In the Ewald summation method there is no interface, 
and the mean potential in the cell is zero. Thus, use of 
Ewald potential to compute charging free energies natu- 
rally conforms to computation of the absolute hydration 
free energy. 

In a cluster, however, the reference potential is not 
zero. Simulations of the vapor-water interface by Sokhan 
and Tildesley[58| suggests that the potential in bulk 
SPC/E water is -12.7 kcal/mole-e relative to the vapor 
value of (zero). Adding this value to our SPC/E esti- 
mate we obtain -256.7 kcal/mole for the hydration free 
energy of the proton in the presence of the potential of 
the phase. This greatly improves the agreement with 
the value of the value of -265.9 kcal/mole obtained by 
Coe and coworkers [HESlEi and -264.3 kcal/mole ob- 
tained by Klots [34| . Consistent with this but probably 
more significantly, the consensus results of TABLE [D] 
suggest a value that is negative and somewhat greater 
than 10 kcal/mole-e in magnitude. 

A classic alternative to the ion-cluster experimental 
studies for obtaining the absolute hydration free ener- 
gies of ions is the tetraphenyl arsonium tetraphenyl bo- 
rate (TATB) hypothesis jSOj. The molecular intuition is 
that these oppositely charged ions have the same non- 
ionic interactions with any solvent molecules. If it were 
precisely true that the absolute hydration free energy 
of these two solutes were precisely the same, then the 
Donnan-like effects discussed above would imply Acj) = 
at equilibrium of this salt between two coexisting flu- 
ids. If true, this would be a satisfactory solution of the 
present problem. But because the values obtained that 
way would not include a contribution from the potential 
of the phase, those values would be different from val- 
ues extracted from the ion-cluster studies analyzed here. 
The computational testing of that TATB hypothesis has 
lead to energetic uncertainties of nearly the same size as 
the energetic contributions of issue here |60|. Therefore, 
the TATB hypothesis must be currently viewed as not 
satisfactorily proven. One important consideration is the 
following: because of the particular electrical asymme- 
try of water molecules, positively and negatively charged 
ions with exactly the same non-ionic interactions with 
water molecules will not have the same absolute hydra- 
tion free energies. It is commonly observed in simulation 



that anions are better hydrated than cations of the same 
non-ionic interactions |9l l6lll62l |: an intuitive view of that 
phenomenon is that water molecule hydrogen atoms, car- 
rying some positive charge, are able to approach anionic 
solutes closely and this is different for the corresponding 
cationic solutes. 

Our computed hydration free energy estimates, espe- 
cially for ion pairs, is in good agreement with experi- 
ments. This agreement should not obscure the approxi- 
mations made in applying a rigorous theory to practical 
calculations. These calculations might be improved in 
several respects. For example, in calculating the excess 
free energy of the cluster and the water ligand, we ne- 
glected packing effects and dispersion interactions. The 
tacit assumption in the analysis is that in composing K n 
(eq. , errors in the cluster and the n water ligands bal- 
ance out. For ionic species this appears to be a good 
assumption, but it does fail for hydration of non-polar 
solutes [63^. We foresee refining these aspects of the cal- 
culation in future studies. 

The neglect of packing contributions is probably more 
easily remedied; see the discussion 2j. Packing contribu- 
tions are expected to be only a few percent of the chem- 
ical and electrostatic effects included already. But those 
packing contributions are expected to be positive. Those 
contributions would probably improve the agreement for 
LiOH and NaOH in TABLE HTH On the otherhand, sev- 
eral additional effects of comparable size would have to 
be included at the same time. For example, the effects of 
anharmonicity on the thermal motion of the clusters is a 
comparable worry. 



VII. CONCLUSIONS 

We find that the absolute hydration free energy, 
/I H+ (w), is somewhat smaller than the acknowledged 
—251 to —264 kcal/mole span of numbers. The present 
analysis suggests that the more negative of these values 
include the potential of the phase, a contribution distinct 
from the absolute hydration free energy. 

Calculations on ion-water clusters can be compared 
properly with experiments on ion-water clusters. Those 
results and comparisons can test the adequacy of such 
calculations. But they don't exclusively test the descrip- 
tion of hydration of the ion; part of those results derives 
from the surface structure of the cluster and that contri- 
bution does not vanish for large clusters. 

This study reinforces the idea that strenuous calcula- 
tion of the absolute hydration free energy, ~p M+ (w), for 
a favorable case such as M=Li could probably reduce 
the computational errors to less than the present exper- 
imental uncertainties. This would permit a more precise 
alignment of available tables of single ion absolute hydra- 
tion free energies. 
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